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In this paper, we proposed a lattice dynamic treatment for the total potential energy for single- 
walled carbon nanotubes (SWCNT's) which is, apart from a parameter for the non-linear effects, 
extracted from the vibrational energy of the planar graphene sheet. The energetics, elasticity and 
lattice dynamics are treated in terms of the same set of force constants independent of tube struc- 
tures. Based upon the proposal, we investigated systematically the relaxed lattice configuration 
for narrow SWCNT's, the strain energy, the Young's modulus and Poisson ratio, and the lattice 
vibrational properties respected to the relaxed equilibrium tubule structure. Our calculated results 
for various physical quantities are nicely in consistency with existing experimental measurements. 
Particularly, we verified that the relaxation effect brings the bond length longer and the frequencies 
of various optical vibrational modes softer; Our calculation provides an evidence that the Young's 
modulus of armchair tube exceeds that of the planar graphene sheet, and the large diameter limits 
of the Young's modulus and Poisson ratio are in agreement with the experimental values of the 
graphite; The calculated radial breathing modes for the ultra narrow tubes with diameter range 
between 0.2 - 0.5 nm coincide the experimental results and the existing ab initio calculations with 
satisfaction; For narrow tubes of diameter 2 nm, the calculated frequencies of optical modes in 
tubule tangential plane as well as those of radial breathing modes are also in good agreement with 
the experimental measurement. In addition, our calculation shows that various physical quanti- 
ties of relaxed SWCNT's can actually be expanded in terms of the chiral angle defined for the 
correspondent ideal SWCNT's. 
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I. INTRODUCTION 

After the discovery of the carbon nanotubes (CNT's) 
in 1991 there have been several methods to prepare 
CNT's with ultra small diameters. Previously, CNT's 
were synthesized in a free spaced 3 - and the quotient of 
narrow CNT's is quite low. Afterwards Z. K. Tang 
et alA&&£ initiated to grew single-walled carbon nan- 
otubes (SWCNT's) inside a channel of zeolite templates. 
The SWCNT's prepared in this way have the diameter 
as small as 0.42 ± 0.02 nm. Actually there are three 
kinds of tubes with different chiralities in this diame- 
ter range: (5,0), (4,2), and (3,3). A fabrication pro- 
cess was further developed to prepare a type of samples 
with (5, 0) and (3, 3) only 8 with a very narrow diam- 
eter distribution. More recently, a stable CNT's with 
diameter 0.3 nm is found inside a multi-walled carbon 
nanotube^ As shown by the density functional studies 9 
that this CNT might be interpreted as the armchair CNT 
(2, 2) with a radial breathing mode (RBM) at 787 cm -1 . 
All these technological improvements stimulated progres- 
sively experimental and theoretical studies on the ultra 
narrow SWCNT's^Ii^ 3 .^^^ 

Since the narrow CNT's possess biggish curvatures, the 
equilibrium geometries would deviate from the ideal ge- 
ometry, i.e., deduced from a seamlessly rolling up pla- 
nar graphitic lattice sheet referring to a chiral vector 



R = niSi + 71,202. 18 Various methods have been de- 
veloped to determine the equilibrium geometry of free 
SWCNT's^iffia Among others, in Ref.M, the au- 
thor calculated the total energy for narrow SWCNT's 
with a first-principles, all-electron, self-consistent local- 
density functional band-structure method. The calcu- 
lated total energy can be parameterized in terms of five 
profitable parameters -i.e., generalized motif variables 
(GMV's)-with the line group symmetry intimately re- 
tained. The equilibrium geometry is obtained by an op- 
timization procedure with respect to these GMV's. In 
Ref. [2(| a molecular mechanical model for the effective 
total potential energy is proposed to calculate the equilib- 
rium structure and the strain energy of achiral SWCNT's 
which is a quadratic form of lattice site displacements 
with respect to the planar graphite sheet. Apart from 
deviations of three bond lengths and three bond angles, 
a pyramidalization angle is introduced to describe the 
energy associated with the curvature. The calculated re- 
sults are consistent with existing numerical results based 
on ab initio calculations. 

Moreover, based upon the relaxed equilibrium geom- 
etry provided by ab initio method^ Milosevic et al^ 
applied the line-group symmetry-based force constant 
method to calculate the vibrational modes of ultra- 
narrow SWCNT's, in which the fitted force constants of 
narrow SWCNT's have been adjusted with respect to the 
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provided cylindrical web geometry. 

Besides, the elastic properties for carbon nanotubes 
have stimulated a large number of experimental and 
theoretical studies^ In particular, the Young's modu- 
lus is a significant physical quantity which reflects the 
striction of equilibrium structure in response to the im- 
pressed load. In Ref. [U, the axis Young's modulus of 
the SWCNT is measured as 1.2 Tpa by atomic force 
microscope technique. By observing their freestanding 
room-temperature vibrations in a transmission electron 
microscope^ the Young's modulus for 27 nanotubes in 
the diameter range 1.0 - 1.5 nm were measured and the 
average Young's modulus is 1.25-0.35/+0.45 Tpa^ An 
interesting fact was acknowledged that the Young's mod- 
ulus of certain nanotubes is even higher than that of bulk 
graphite. 

There have been a large number of theoretical cal- 
culations for the Young's modulus and the Poisson 
ratio. The ab initio method , 21 i 27 i 28 tight-binding 
approximation ) 29 i 30 i 31 and force constant model 32 gave 
the results of Young's modules with a range from 1 to 
1.24 Tpa, and Poisson ratio in a range 0.16 - 0.27. It fur- 
ther shows that the calculated results vary slig htly w ith 
different radii and chiralities. However, in Refs 1 3 3U 3^f35l . 
[H, the calculated Young's modulus or the Poisson ra- 
tio are found to be chirality and tube radius-dependent. 
In Ref HH, the Poisson ratio of armchair tubes has the 
value larger than that of planar graphene sheet while for 
zigzag tubes smaller, and its planar sheet limit is 0.21. 
In Refs [34] and HH, the Poisson ratio decreases with in- 
creasing the tube diameter with the large diameter limit 
as 0.16. Anyhow, to our knowledge the resolution of 
the experimental measurements for the Young's modu- 
lus and Poisson ratio is still not enough satisfactory yet 
and the varies correspondent theoretical calculations are 
even more diverse. 

As we have seen in above, to figure out the total po- 
tential energy of SWCNT's as the functional of cylindri- 
cal lattice configuration is the key issue underlying the 
investigation of the equilibrium structure, strain energy, 
Young's modulus and lattice dynamics of ultra-narrow 
tubes. 

In this paper, we propose that the total potential en- 
ergy for SWCNT's can be extracted from the vibrational 
energy of the planar graphene sheet with five force con- 
stants. In which the bi-linear forms of the lattice site 
displacements are accounted with the equilibrium lat- 
tice configuration of the planar graphite as the reference 
point, while one of the potential term for the twist mo- 
tion needs to be improved by including a quartic term to 
account the non- linear effect. One of the advantages of 
our proposal is that the equilibrium structure, the elas- 
ticity, and the lattice dynamic properties for all kinds of 
SWCNT's are calculated in terms of the same set of six 
(5+1) force constants, which essentially reflects an uni- 
fied microscopic mechanism founded on the five forms 
of motion of the graphite- type cylindrical lattice sheet. 
Based upon the proposal, we investigate systematically 



the relaxed lattice configuration for narrow SWCNT's, 
the strain energy, the Young's modulus and Poisson ra- 
tio, and the lattice vibrational properties respected to 
the relaxed equilibrium tubule structure. The calculated 
results are nicely in consistency with the existing exper- 
imental measurements. 

In particular, our calculation reveals a generic feature 
that the relaxation effect exhibits itself in stretching the 
bond lengths as well as softening the mode frequencies 
for all kinds of SWCNT's; It provides a kind of evidence 
that the Young's modulus of armchair tubes exceeds that 
of the planar graphene sheet, i.e., when the radius in- 
creases, the calculated values of the Young's modulus 
for armchair tubes approach that of the graphite from 
above and those of zigzag tubes from below. Mean- 
while, the large diameter limits of the Young's modulus 
and Poisson ratio are in agreement with the experimen- 
tal values of the graphite; The calculated RBM for the 
ultra narrow tubes with diameter range between 0.2 - 
0.5 nm coincide the experimental results^ and the exist- 
ing ab initio calculations^ with satisfaction. For narrow 
tubes of diameter 2 nm, the calculated frequencies of op- 
tical modes in tubule tangential plane as well as those of 
RBM are also in good agreement with the experimental 
measurement^ The honeycomb lattice symmetry-based 
chiral expansion for physical quantities of SWCNT's with 
ideal geometry 3 ^ can be retained for those of relaxed 
SWCNT's by changing the chiral angle into its counter- 
part defined in tubes of ideal geometry; As expected, for 
tubes with diameter larger than approximately 0.8 nm 
where the relaxation effect dies away, the calculated fre- 
quencies and the sound velocities for various vibration 
modes are in accord with those of the tubes of ideal 
geometry^ 

In Section II, the total potential energy for SWCNT's 
is presented, and the relaxed equilibrium geometry and 
correspondent strain energy are discussed. Section III 
is devoted to the lattice dynamics for ultra-narrow 
SWCNT's, while Section IV for the Young's modulus and 
Poisson ratio. Finally, the conclusion is in Section V. 



II. THE RELAXED EQUILIBRIUM 
STRUCTURE 

A. The geometry 

We recall the geometrical description for the relaxed 
equilibrium structure of SWCNT's particularly for the 
ultra-narrow tubes following Ref. ■10). As for the cylindri- 
cal coordinate description shown in Fig.[TJ we first set the 
z axis along the tube axis, while the x axis is fixed pass- 
ing through one of the A atoms. We further introduce 
a planar image of the cylindrical lattice sheet by unfold- 
ing the cylindrical tube into a planar lattice sheet, which 
is actually a tangential plane with the tangent line as a 
generatrix passing through the A atom, and set atom A 
as the origin of the plane with coordinates as (r^,0,0), 
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r s = r, 4> = ma + ^27rr + a', z = to/i + /i' (for B). 

While the corresponding coordinates in the image plane 
have the form as 



f imag0 (m,/,s) = {r s ,cjrr s ,z) 



(2) 



with the ^-components of both A and B atoms as con- 
stants respectively. 

Therefore, the six GMV's r, ta, a', h', a, and h provide 
a complete description for the geometry of the lattice 
structure of SWCNT's, in which all the A and B atoms 
are allowed to sit on different cylindrical surfaces with 
radii ta and r respectively. 



FIG. 1: Sketch of carbon atoms on a cylindric surface. 



B. The total potential energy 



where ta is the radius of atom A. The planar lattice 
sheet is constituted of parallelogram primitive cells with 
di and 3,2 as its basic vectors. The deviations of CNT's 
with relaxed geometry from that of ideal geometry can 
be illustrated most conveniently by the primitive cells 
on the image plane. For that of ideal geometry, the basic 
vectors are of equal magnitude with an intra-angle of 7r/3 
between themselves, while for that of the relaxed geom- 
etry, their magnitudes are no more equal and the angle 
can be deviated from tt/3. Meanwhile, the location of the 
in-cell B atom relative to A atom could also be deviated 
from that of the ideal geometry. 

Now we introduce the chiral vector R = n±ai + nidi 
initiated from the atom A along the y axis in the im- 
age sheet, and R is the image circumference of cylin- 
drical lattice sheet. We then introduce a screw vec- 
tor H = p\d\ + pidi with integers p\ and pi satisfying 
p\n% — pin\ = N and N is the largest common divisor of 
the integer pair (ni,ni). Exactly following RefsllfJ and 

[Lsl . H and constitute the basic ingredients of the gen- 
eralized motif cell, upon which a generalized screw oper- 
ation S(a, h) and a generalized AT-fold rotation operation 
Cjv can be established. The former is a rotation a around 
the tube axis with a simultaneous translation h along the 
tube axis, the latter a rotation of -3 around the tube axis. 
A generalized motif cell (to, I) can be generated from the 
(0, 0) cell by successive operations of S(a, h) and Cn for 
to and I times respectively. For the in-cell geometry of 
a generalized motif cell, originated from atom A the po- 
sition of atom B can be fixed by a rotation a' around 
the tube axis and two simultaneous translations h! and 
(r — ta) along the tube axis and the tube radius respec- 
tively. 

As a result, the coordinates on the cylindrical lattice 
sheet for each site (m, I, s) with s = A,B can be ex- 
pressed as 



r(m, I, s) = (r s cos <f>, r s sin 0, z) 



(1) 



In fact, in sense of quadratic approximation the vi- 
bration potential of the single layer graphite lattice sheet 
provides a proper description of various modalities of lat- 
tice motion. We, therefore, propose that the total po- 
tential energy of SWCNT's can be extracted from the 
vibration potential of graphene lattice sheets as follows. 
Firstly, the total potential energy contributed by carbon 
atom I, i.e., atom (0,0,^4), is composed of the following 
five terms in which a local planar graphene sheet is taken 
as the reference point, 



U A = Vl + V 8 i + V BB + V r , 



Vu 



(3) 



where Vi, V s i are the potentials of the spring forces 
between the nearest-neighbor and the next-nearest- 
neighbor atom pair respectively as shown in Fig.[TJ Vb b is 
the energy associated with bond angle variations, V rc de- 
scribes out-of-surface bond bending, i.e., a kind of strain 
force on atom i in out-of-surface direction by its three 
nearest neighbors in curling process, and Vt w is the twist 
potential energy. 



ta , 4> = met 



A 



2tt, z = mh (for A), 



Vi = 

V sl = 

V BB = 

V rc = 

v tw = 

x 1 - = 

m = 



[(Af - - An) • e{ 



4 ^ 

i=2 



, 10 

3=5 



I 12 



?; 12 



2 u 



(COS &r) — COS 0q) 2 



J2 - 3Ar! I • e 1 



tE/(^). 

[(An - Afj) - (Arv - Ar r )] ■ e{ k , 



(4) 
(5) 
(6) 

(7) 
(8) 



K ah x 4 , if Id < 



2K a 



IK„ 



4 



In Eqs (gl) to © i = 2, . . . , 4 and j = 5, . . . , 10 are 
the nearest and next-nearest neighbors of atom 1 respec- 
tively, Qrj with rj = 1,2,3 represent the three bond angles 
with atom 1 as the common apex while 0o = ^ is the 
bond angle in the graphic plane, represents a pair 

of atoms nearest neighbored to atom 1 with k being the 
third of its nearest neighbors and the pair is the 

image of referring to a C2 rotation around the axis 
in e[ k . Moreover, 



AFj =fi — Ft 



(9) 



is the displacement of the carbon atom i from the planar 
graphene sheet r, to its counterpart in the tubal lat- 
tice sheet Fi, in which i stands the site index (m,l,s) 
defined on the cylindrical lattice sheet. In addition, 

e ls = r^r^rj is the unit vector pointed from atom 1 to 

atom i in the graphite sheet, e[ and e{ k are vertical unit 
vectors of the planar graphene sheet located at the site 
% and the middle point of sites 1 and k of the graphene 
sheet respectively. All the unit vectors are defined on the 
local reference graphene sheet and are introduced to keep 
the rigid rotational invariancc symmetry. The above five 
potential terms in principle cover the main features of 
the modalities of lattice deformation suggested by Lianxi 
Shen et a/. ) 20 ' 34 Chang et al^ in their mechanical model. 

Correspondingly, that part of total potential energy in 
association with carbon atom B (0, 0, B) has a similar 
expression as that of Ua, but with its geometrical pa- 
rameters replaced by the counterpart of atom B. 

Secondly, we express the lattice displacements in Ua 
and Ub in terms of the generalized motif variables. It 
spells that the locally constructed (0, 0) generalized motif 
cell with a pair of carbon atoms A and B can be contin- 
ued to form a whole seamlessly closed cylindrical lattice 
sheet in which the generalized screw as well as rotation 
symmetries are inherited. As a result the total potential 
energy of SWCNT's has the form as a simple multiple of 
those associated with a pair of nearest neighboring atoms 
A and B, i.e., Ua and Ub respectively, 



U = JV*(U A + U B ) 



(10) 



where M is the number of generalized motif cells of 
SWCNT's in consideration. The redundant counting in 
the total potential energy contributed by different car- 
bon atoms can be accounted by a proper resetting for 
the values of the relevant force constants. We further 
address that this sort of parameterization procedure can 
be carried out for all kinds of SWCNT's including both 
chiral and achiral tubes. 

The values of force constants are fc; = 305 N m _1 , 
k sl = 68.25 N m -1 , k BB = 1-38 x lO" 11 erg, k rc = 14.8 N 
m -1 , k tw = 6.24 N m -1 , and K ah = 2.5. The frontal 
five force constants come from vibration energy of planar 
graphene sheet 39 . Since the twisting deformations in- 
crease substantially with decreasing the tube diameters, 
for example the twisting angle for the C-C bond can take 
the value as large as 26° for the narrow SWCNT (6,3), 



TABLE I: The calculating results of the five GMV's with 
fA = r and strain energy for narrow SWCNT's. 



(m , U2) 


r (Aj 


OL (TcLQJ 




cm. (rsLCi) 


n [A) 


i2j s \e v ) 




1. f 4o 


1. J.UOO 


u.uuu 


1.0 / Uo 


l.ZOO 


1 CCQC 

i.ooyo 


a 1 
(A JJ 


1 K1 7 


1 nni q 


u.oo / 


A ^09^ 


n ^sn 

U.OoU 


1 7p:qc 

1. / oyo 


(4, °J 


1 ri77 

i.y/ / 


n 7QC A 

U. f o04 


U.olo 


U. /o04 


nQQ 
z.Uoo 


l.ODOD 


V ' z ) 


niQ 
z.Ulo 


U.o4dU 


ft 1 AS2 


9 A 111 


n /ioi 
U.4yl 


n 7077 
U. 1 y 1 t 




9 1 9Q 

z.izy 


n 7^1 9 

U. 1 01Z 


U.41D 


A 0*31 7 


n a 


n si 71 


(■5, 6) 


9 981 
Z.Zol 


U. ( Uoo 


n nnn 
u.uuu 


1 fi/1 79 
1.U4 / Z 


1 9Q2 
l.ZOO 


n A Pi99 
U.40zZ 


(5, 0) 


2.283 


0.6283 


0.652 


0.6283 


2.069 


0.6484 


(4, 2) 


2.333 


0.6788 


0.250 


2.0127 


0.803 


0.4781 
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0.6075 
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0.376 


0.4606 
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0.6000 
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1.7827 
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0.3083 


(6, 0) 
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0.5236 
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0.3759 


(5, 2) 


2.673 


0.5652 
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2.6560 


0.339 


0.3163 



TABLE II: Contributions of five potential energy terms to the 
strain energy for tube (3, 3). 



(3, 3) 



v s . 



Vbi 



V r , 



v tu 



(eV) 0.02934 0.07191 0.02972 0.30179 0.01948 



an additional force constant for the anharmonic improve- 
ment as K a h in Eq. ([5]) has to be introduced. We notice 
here that the differences in lattice dynamics between such 
two types of systems are essentially due to their different 
spatial geometry, which are characterized by their own 
equilibrium lattice structure. 



The relaxed structure of SWCNT's 



We introduce the total potential energy per atom as 



E, 



E = E{r, r A ,a, h, a', ti) = (U A + U B )/2 



(11) 



By minimizing E(r,rA, a, h, a 1 , h') with respect to the 
six GMV's, we obtain the optimized equilibrium geome- 
try described by GMV's (f, r~A, a, h, a', h') as well as the 
strain energy E s = E(r,r~A,a,h,a',h'). It is interesting 
to address that the optimized equilibrium lattice config- 
urations for all tubes are shown to be of t~a — f. This is 
because that the total energy U keeps unchanged when 
the tube is rotated upside down with atoms A and B mu- 
tually permuted, so that the optimized geometry must 
exhibit the C2 symmetry. The reflect symmetry ah with 
respect to the cross section is also kept for achiral tubes. 

The calculated various lattice structure properties and 
the strain energies for narrow SWCNT's with diameter in 
[0.26, 0.5] nm are listed in Table HI where the calculated 
GMV's are in good agreement with those by ab initio 
calculations. 10 Moreover, it is shown that the curvature 
energy associated with the fourth potential energy term 
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FIG. 2: The strain energy as the function of radius for 
all tubes in [1.3, 10.0] A. It is compared with the ab initio 
calculations . ' 21 



in Eq. ([7]) dominates the total strain energy, while the 
corresponding contributions from the bond variation are 
usually much less than those from the bond angle varia- 
tions (see Table HI]) . 

In Ref. HH, it is reported that the strain energy various 
with tubule radius as -K by applying the Tersoff-Brenner 
interatomic potentials. We here plot the strain energy 
versus radius in Fig. [2] which is comparable with the ex- 
isting resul t 21 ' 36 and can be well fitted by E s = % + 
in consistency with the result of Ref. [2(| The coefficient 
of ^-i.e., C2~is a bit large. That is because, in Ref. I20I 
there is no such twisting energy term as we have. 

Our calculation shows that the radius and C-C bond 
lengths of CNT's of the optimized lattice configuration 
are always larger than those of the correspondent CNT's 
of the ideal geometry. We representatively display the 
three bond lengths in unit of the graphene bond length 
0.142 nm for (2n,n) tubes with the relaxed geometry as 
well as ideal geometry in Fig.[3^i. The relaxation, resulted 
from the optimization procedure, would disentangle the 
warping-induced tension and make the inter atomic bond 
length getting longer. To be precise, the bond lengths of 
the optimized lattice configuration decrease with increas- 
ing diameter and keep their values always larger than 
that of the graphene sheet. In contrary, for tubes of 
ideal geometry the bond lengths increase with increas- 
ing diameter and keep their values always smaller than 
that of the graphene sheet. As we can see in Fig. [3k,, 
one of the bond lengths of the tube (4, 2) is 0.139 nm for 
ideal geometry while 0.157 nm for case of relaxed geom- 
etry which is about 11% more than that of the graphene 
sheet. 

As shown in Fig. [5b, deviations of the bond angles 
from that of the graphene are really considerable for nar- 
row tubes with diameter smaller than 0.8 nm. On the 
other hand, as can be seen in the same figure, this effect 
decreases quickly with increasing the tube diameter. 
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FIG. 3: C-C bond lengths reduced by that of graphite and 
bond angles for different chiral tubes (2n, n) with the relaxed 
geometry and ideal geometry. 



III. CALCULATION RESULTS AND 
DISCUSSION OF PHONON DISPERSION 

In Sec. II, the relaxed equilibrium lattice configura- 
tions of SWCNT's are derived by an optimization proce- 
dure applied to the proposed total potential energy. The 
vibrational potential of SWCNT's is actually caused by 
the deviation of the lattice sites with respect to the re- 
laxed equilibrium lattice configuration. It is obvious that 
the vibrational potential of SWCNT's would take the 
same form as the total potential energy, whereas the equi- 
librium position of the lattice site in the planar graphene 
sheet is therefore replaced by those in the relaxed cylin- 
drical lattice sheet, meanwhile, the non-linear term is no 
more needed. The detailed modifications are as follows. 

(1) Afi — Ti — rto is renamed by convention as the 
vibrational displacement Ui with rio the local planar 
graphene sheet to be replaced by the relaxed equilibrium 
position of atom i in SWCNT's; 

(2) cos Oo in Eq. ([6]) is no more the bond angle of 
graphene sheet as |7r but the correspondent bond angle 
of the relaxed equilibrium configuration which depends 
on distinct adjacent atom pairs nearest neighbored to 
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the common apex atom i; 

(3) and e/- are kept to be the unit vectors pointed 
from atom i and j to atom 1 respectively. e{ is replaced 

by e[ c — — ~;' 4 = 2 — and eVl becomes now the unit vec- 

lL, 4 =2 fio l < 

tor along the radial direction of the middle point of atoms 
1 and fc. We stress that all these unit vectors are defined 
on the cylindrical lattice sheet with optimized equilib- 
rium geometry, which taken a crucial role to keep the 
rigid rotational invariance for the vibrational potential 
of relaxed SWCNT's; 

(4) The coefficient K a h is taken to be zero. 

We emphase here that the five force constants of the 
above proposed quadratic vibrational energy are the same 
as those in Sec IIB. They are applied to all kinds of 
SWCNT including those with small radii. This is be- 
cause that the rigid rotation symmetry of the SWCNTs 
is precisely kept. And the curvature effect for the vibra- 
tional modes has been properly taken care by the above 
introduced quadratic expression. 

It is straightforward that such obtained vibrational po- 
tential can be again parameterized in terms of GMV's. 
Taking the advantage of GMV description for the cylin- 
drical lattice configuration, the underlined generalized 
screw and rotation symmetries ensure that the lattice dy- 
namic equation even for tubes with relaxed geometry can 
be decomposed into a six-dimensional eigenvalue problem 
and becomes treatable. 

We calculate, for tubes with diameter from 0.28 nm to 
2.5 nm, all the Raman- and Infrared-active modes^ and 
the velocities for the twisting modes (TW) and longitude 
acoustic modes (LA). 
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FIG. 4: The radius dependence of relative frequency change 
due to the relaxation effect for the optical mode along the z 
axis. The small waviness of the line shows the dependency 
of the frequencies on the chiral angle. For chiral and zigzag 
tubes this mode is the Raman-active mode. 

The relaxation effect for the non-zero vibrational 
modes of the tubes with relaxed geometry exhibits itself 
as the frequency lowering (mode softening) in compari- 
son with their counterpart of the tubes of ideal geometry. 



TABLE III: The values of two optical modes in the tangential 
plane from our theory and compared with the experimental 
results.— 



Optical mode 


(11, 10) 


(17, 9) 


(27, 4) 


(15, 14) 


Our theory 


1597.3 


1600.5 


1602.5 


1602.1 




1600.7 


1603.6 


1605.5 


1603.7 


Experiment 37 


1566(7) 


1572(13) 


1577.5(10) 


1573(7) 




1593.5(6) 


1591(8) 


1593(5) 


1592(5) 


TABLE IV: The values of RBM from 


our theory and com- 


pared with the experiments or other calculations. 




RBM 


(2, 2) 


(5, 0) 


(3, 3) 


(4, 2) 


Our theory 


804.3 


561.3 


550.8 


538.6 


comparison 


(787) a 


(550) b 


/ 


(510) 6 


RBM 




(11, 10) 


(16, 7) 


(15, 6) 


Our theory 




160.0 


142.5 


155.2 


comparison 


(169.5(7)) c 


(154(5)) c 


(166(l)r 



a The ab initio calculations from Ref. @; 6 The experiment 
results from Ref. 3 c The experiment results from Refs [3?] 
and lU 



The effect of the relaxation on the optical mode along 
the z axis is shown in Fig. [4j The relaxation softens this 
mode considerably for narrow tubes. For example for the 
tubes with radius as narrow as 2 A, this mode is softened 
as much as 8% in reference to that of the ideal geometry. 

The calculated frequencies for the optical mode in 
the tubal tangential plane coincide excellently with the 
experimental measurements on tubes (11,10), (17,9), 
(27, 4), (15, 14) ,21 see Table EH For the RBM, the calcu- 
lated results from our model arc in good agreement with 




K[2jt] 

FIG. 5: Phonon dispersion curves for SWCNT (4,2) in the 
representation of helical and rotational quantum numbers. 
The flexure mode — i.e., the transversal acoustic mode — shows 
up at (k, n) — (a, 1), where a is in the unit of 2ir. It is double 
degenerate. 
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TABLE V: Polarization vectors u = (u(A),u(B)) as functions of r (r £ [1.3, 10.0] A) and 0. Where u{A) and u(B) indicate the 
displacement vectors of atoms A and B in the (0,0) unit cell respectively. For the three modes in this table, u r (B) — u r (A), 
u<f,(B) = -u^{A), u z {B) = -u z {A). 



(k,ti) = (0,0) 




Vector(0) 


Mr) 






Ri (e r AC) 


u r (A) 


fair) 


fair) 


- 0.7072 °'T 6 






U4>{A) 


/i (r) sin 30 


Mr) 


0.0212 | 0.1111 
— ~r 5 — 

r t 






u z (A) 


/i (r) cos 30 


Mr) 


0.0344 | 0.2167 
r r -z 




R 2 (e^ OP) 


U r (A) 


fi (r) sin 38 


Mr) 


0.0211 0.1064 
— " ' — 2 — 






u^A) 


Mr) +/i(r) cos 12(9 


Mr) 


= 0.6975 + 2%2, / x ( r ) 


- 0.0084 M2« 




u z {A) 


/i (r) sin 60 + / 2 (r) sin 120 


Mr) 


= 0.1642- ai^a, /a(r ) 


= 0.0283 + 2il|12 


R 3 {e z OP) 


Ut(A) 


/i (r) cos 30 


Mr) 


_ 0.0348 0.2189 






u^{A) 


/i(r) sin 60 + / 2 (r)sinl20 


Mr) 


- 0.1643 Mpa,/ 2 ( r ) 


- 0.0283 + 2il|10 




u z {A) 


/o(r)+/i(r) cos 120 


Mr) 


= 0.6975 + 2^S£ 1 / 1 (r-) 


- 0.0084 2^M£ 



the correspondent results from the experiments on tubes 
(5,0), (4,2), (11,10), (16,7), (15,6) listed in Table [TV] 
Also for tube (2, 2) the calculate result is consistent with 
that of the ab initio calculation 9 (see Table IIV[) . It can 
be seen moreover from Table IIVI that the smaller tube 
has the larger RBM frequency. This is consistent with 
some other calculation resultsj&i 1 - 

We take SWCNT (4, 2) as a typical example of the 
narrow tubes. Its dispersion is shown in Fig. [5] with the 
helical and rotational quantum numbers k and n as n = 
0, 1 and n S [0, 7r]. In the neighborhood of (k, n) — (0, 0), 
there are two zero-frequency modes, TW and LA, with 
acoustic velocities as Ctw = 13.6 km s -1 , and Cla — 
20.7 km s _1 respectively. It is interesting to find that the 
Ctw is much smaller than that in the graphite plane. 
This effect is expected to be measured in experiments 
with technical developments. The frequencies of RBM, 
e r optical (OP) mode, OP mode and e z OP mode 
are calculated as 538.6, 1081.5, 1544.2 and 1487.7 cm" 1 
respectively. At (n,n) — ±(a, 1) with a = 2.01 radian, 
there are two transverse acoustic modes- i.e., the flexure 
modes- with the parabolic dispersions as lu 2 — /3 2 (K=Fa) 4 
in the low-frequency limits. They are the consequence of 
the rotational invariance of our vibration potentia h 38 ' 42 

For tubes with diameter larger than approximately 
0.8 nm where the curvature effect dies away, the calcu- 
lated frequencies (Fig. |4|) , the sound velocities for various 
vibration modes, and even bond length (Fig. [3]) etc. are 
in accord with those of the tubes of ideal geometry^ We 
notice that the slight difference in force constants be- 
tween the present paper and Ref. [38| take no substantial 
effect on the lattice vibrational properties. 

Although the lattice structure for the relaxed nar- 
row SWCNT's deviates from that of the ideal geome- 
try considerably, the chiral index (ni, n-i) for the relaxed 
SWCNT's retains the symmetry information inherited 
from the hexagonal planar lattice sheet. As one might ex- 
pect, our calculation verifies that various physical quan- 
tities can be expanded in terms of cos3n# and sin3n#, 
n = 1,2,..., with the same expansion forms as those 



of the ideal geometry, in which 9 is the chiral angle de- 
fined for the correspondent tube with ideal geometry as 
9 = arctan „^ n2 . As an example, we show our fitted 
polarization vectors of RBM and two optical modes in 
tangential plane at (K,n) — (0,0) in Table fVl The corre- 
spondent coefficients are shown almost the same as those 
in Table IV of Ref. M. 



IV. THE YOUNG'S MODULUS AND THE 
POISSON RATIO 

The relaxed lattice configuration provides itself as the 
proper minimum of the total potential energy. In partic- 
ular the equilibrium lattice configuration for the loaded 
SWCNT's can also be performed, which results that the 
calculation for various strain-induced physical quanti- 
ties such as Young's modulus and Poisson ratio can be 
straightforwardly carried out. Among the GMV's, h and 
h' are the variables of the longitudinal degrees of free- 
dom, while r, a and a' are variables of peripheral de- 
grees of freedom. However, the experimentally measured 
longitudinal stretching should correspond to the variable 
h while the variable h! essentially describes the micro- 
scopic in-cell displacement. The Young's modulus for 
the SWCNT's are introduced as: 



11 zero stress except the (z,z) component 

where e±i = ^ being the longitudinal strain, and the 
derivatives are taken with all the other GMV's r, a, a' 
and hi being freely relaxed which corresponds to the 
boundary condition that all the components of stress 
tensor are zero except the (z, z) component. We stress 
that such a calculation of the Young's modulus for the 
cylindrical lattice sheet is consistent with the classical 
definition of Young's modulus in the continuum limit 42. 
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TABLE VI: The large diameter limits of Young's modulus and 
Poisson ratio from our calculation and other existing results 
are compared with the experiment data. 





Experiment 4 * 


Our work Theory" Theory 6 


Theory c 


Y (eV) 


56.43 


56.3 60 69 


53.3 




0.17 


0.166 0.19 0.25 


0.277 



a The ab initio method^ 6 The tight-binding 
approximation; 2 - c The force constant model^ 2 - 



Accordingly the Poisson ratio 



£22 
£11 



zero stress except the (z, z) component 



(13) 



with £22 being the peripheral strain induced under the 
same condition as explained above. 

The calculated Young's modulus and Poisson ratio for 
various narrow SWCNT's can be well fitted as the func- 
tions of radius r and chiral angle 9, 



Y 



56.3 



15.1 + 21.8 cos 69 



(14) 



0.235 + 0.246 cos 66» 
H = 0.166+ 5 , (15) 



where Y is in the unit of eV, r is in the range of 
[1.3,10] A, 9 is the chiral angle defined for the corre- 
spondent SWCNT with ideal geometry as explained in 
the previous section, and the relative fitting error is kept 
less than 1 x 10~ 3 . As shown in Eqs (|T3)) and (TPS")) , 
the Young's modulus and Poisson ratio in large diameter 
limit have the values of 56.3 eV and 0.166 respectively, 
which are in good agreement with the experimental re- 
sults of graphite^ 4 - (see Table IVI[) . In the case of narrow 
tubes, they are significantly chirality-dependent. For ex- 
ample, the Young's modulus of three SWCNT's (3,3), 
(4,2), and (5,0) calculated as 57.74, 54.94, and 49.26 
are considerably different from each other, although they 
have almost the same diameter. 

Furthermore, we display the radius dependence of the 
Young's modulus in Fig. [6] for two species SWCNT's: 
armchair tubes (n, n), and zigzag tubes (n,0). It shows 
clearly that the Young's modulus of armchair tubes de- 
creases with increasing of tube diameter and always takes 
the values larger than that of planar graphite. While for 
zigzag tubes it exhibits a contrary behavior, i.e., increases 
with increasing diameter and keeps its value always below 
that of graphene sheet. It was reported in Ref. [26j that 
among 27 samples of carbon nanotubes, some of them are 
measured with the Young's modulus higher than that of 
bulk graphite. While, to our knowledge, most of the ex- 
iting theoretical estimations^ can only provide values of 
Young's modulus less than that of graphene sheet for any 
kind of SWCNT's. Our calculation provides an evidence 
that armchair tubes are stiffer while the zigzag softer due 
to the chirality-dependence. 



T3 

o 
E 



60 



56 



52 



48 



44 



40 



Young's modulus 



(n, n) Energy method 
(n, n) Velocity method 
(n, 0) Energy method 
(n, 0) Velocity method 



4 6 
radius [angstrom] 
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10 



FIG. 6: The Young's modulus for armchair and Zigzag tubes. 
The lines are calculated by the second derivatives of energy 
per atom with respect to the uniform strain along the tubule 
axis direction. And the dashed lines are calculated by Eq (|16|) 
with velocity method. 
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o 

in 
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FIG. 7: The Poisson ratio for different armchair and zigzag 
tubes. 



It is known that following the elasticity theory of the 
continuum medium, the Young's modulus as well as 
Poisson ratio can be derived from the relevant sound 
velocities Mdl 



Y 



fi = 0.5 * 



Vl 



1, 



(16) 



(17) 



where the two expressions are also referred to the bound- 
ary condition of only the (z, z) component of stress tensor 
being non-zero. Therefore, we may calculate the Young's 
modulus and Poisson ratio from the relevant sound ve- 
locities obtained in the last section shown in Fig. [6] for 
comparison. We emphasize that the Young's modulus 
(Poisson ratio) calculated from the strain energy follow- 
ing Eqs fT2")) and (fl"3"|) and those calculated from the 
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sound velocities following Eqs (fT6|) and (fT7f are essen- 
tially two distinct independent calculations. The former 
is in sense of a microscopic calculation following from the 
total potential energy Eqs ||3J|-||HJ), while the latter is a 
sort of macroscopic calculation with the sound velocities 
as its input although which is provided by the micro- 
scopic vibrational energy (introduced in section III). We 
can see from Fig. [13 for each of the two species of SWCNT 
(n,n) and (n, 0), the Young's module produced by the 
two kinds of calculations meet to each other closer and 
closer as the radius increases and have the correct r — > oo 
limit as the experimental Young's modulus of graphite. 
Such a result arises from the intrinsic consistency be- 
tween our total potential energy expression and that of 
the vibrational energy. The difference between the two 
kinds of calculations becomes apparent when tube radius 
decreases. We understand such a phenomenon as that 
the calculation following from the theory of elasticity is 
macroscopic in nature which would be failure for ultra 
narrow tubes. 

We show the calculated Poisson ratios for armchair 
{n,n) and zigzag (n, 0) tubes in Fig. [7] as the function of 
radius. The values of the Poisson ratio are in an order 
as that those of zigzag are the largest and the armchair 
the smallest. This is consistent with the correspondent 
result of ab initio calculation^ 



V. CONCLUSIONS 

In this paper, we proposed a lattice dynamic treat- 
ment for the total potential energy for SWCNT's which 



is, apart from a parameter for the non-linear effects, 
extracted from the vibrational energy of the planar 
graphene sheet. Based upon the proposal, we investi- 
gated systematically the relaxed lattice configuration for 
narrow SWCNT's, the strain energy, the Young's mod- 
ulus and Poisson ratio, as well as the lattice vibrational 
properties respected to the relaxed equilibrium tubule 
structure with the same set of force constants indepen- 
dent of tube structures. In particular, with the applica- 
tion of GMV's we not only successfully extend the locally 
introduced-total potential energy to the whole tube, but 
also make the lattice dynamics for SWCNT's with re- 
laxed geometry becomes treatable. Comparing with the 
corresponding results obtained from the tubes with ideal 
geometry, we verified that the relaxation effect brings the 
bond length longer and the frequencies of various optical 
vibrational modes softer. Moreover, our approach not 
only deduces the proper large diameter limit values of 
Young's modulus and Poisson ratio, but also provides a 
strong evidence that SWCNT's with different chirality 
could be stiffer or softer than the planar graphene sheet. 
The calculated RBM and optical modes in tubule tangen- 
tial plane are also in good agreement with the experimen- 
tal measurements. In addition, we realized that the chiral 
symmetry-based general expansion formulae for different 
physical quantities survive for SWCNT's with relaxed ge- 
ometry in the sense to replace the chiral angle by its 
counterpart defined in tubes of ideal geometry. 
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